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ABSTRACT 

Previous ground based observations of the Seyfert 2 galaxy Mrk 78 revealed a 
double set of emission lines, similar to those seen in several active galactic nuclei 
(AGN) from recent surveys. Are the double lines due to two AGN with different 
radial velocities in the same galaxy, or are they due to mass outflows from a single 
AGN? We present a study of the outflowing ionized gas in the resolved narrow- 
line region (NLR) of Mrk 78 using observations from Space Telescope Imaging 
Spectrograph (STIS) and Faint Object Camera (FOC) aboard the Hubble Space 
Telescope(HST) as part of an ongoing project to determine the kinematics and 
geometries of AGN outflows. From the spectroscopic information, we determined 
the fundamental geometry of the outflow via our kinematics modeling program 
by recreating radial velocities to fit those seen in four different STIS slit positions. 
We determined that the double emission lines seen in ground-based spectra are 
due to an asymmetric distribution of outflowing gas in the NLR. By successfully 
fitting a model for a single AGN to Mrk 78, we show that it is possible to explain 
double emission lines with radial velocity offsets seen in AGN similar to Mrk 78 
without requiring dual supermassive black holes. 
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Introduction 



Narrow line regions of Seyfert galaxies contain photoionized gas due to the nonstellar 
continuum emitted by the AGN and are structured rough ly in the shape of a bicone, wher e 
the apex of the bicone is located in the center of the AGN (|Poggelll988l ; ISchmitt et al.lll994j ). 
This likely indicates that the ionizing radiation from the nucleus is being collimated by an 
absor bing material, which could be a torus-like structure (lAntonuccilll993l ) or perhaps a disk 
wind ( lElitzur fc Shlosmanll2006l ). However, the source of the NLR gas is not well understood. 



An AGN of particular interest is Mrk 78. Classified as an SB galaxy in the NASA/IPAC 
Extragalactic Database (NED), it is also a Seyfert 2 galaxy because only narrow (full width 
at half-maximum [FWHM] < 1000 km s _1 ) emission lines are present in i ts optical spectra. 
With a redshift of z — 0.036793 based on Ca T stellar absorption lines ( INelson fc Whittle 
19951 ). Mrk 78 is at a distance of ~ 150 Mpc (for Hq = 73 km s -1 Mpc -1 ); at this distance, 
1" corresponds to a transverse size of ~730 pc. Mrk 78 has been extensively studied due to 
a visual correlation between the its radio source and extended n arrow-line region (ENLR) 
((Whittle et aJbood Iwhittle fc Wilsonlbooi Iwhittle et ajjooi ). HST FOC images of Mrk 
78 in [O III] show that its narrow-line region (NLR) consists of two highly luminous regions 
extending over ~ 5" (~ 3.6 kpc) that is in a biconical geometry also seen in several other 
Seyfert 2 galaxies ( ISchmitt et al.l 120031 ) . Thus, Mrk 78 makes for an interesting addition to 
our ongo ing study to determine the geometries of the outflows and fueling flows in the NLRs 
of AGN foas et al.l 120051 bood : ICrenshaw etaJboiOal iFischer et aDboioh . 



A specifically interesting aspect is that ground-based spect roscopic observations of Mr k 
78 have shown that its narrow emission lines are double peaked (jSargentlll972l ; lAdamslll973l ). 
Ground-based spectra at high spectral resolution depi ct the [O IHl emission-line profile as two 
large , well-defined peaks separated by ~ 800 km s _1 ({Whittle et al.lll988t INelson fc Whittle 
19951 ). As of late, there has been considerable activity around the study of double emission 
lines as they have been used as an indication for the presence of multiple AGN in a single 
galaxy. As such, Mrk 78 would seem an excellent candidate to test for multiple AGN as it 
is nearby, can be studied in detail, and contains dual-peaked lines with a velocity separation 
similar to those seen in other ex amples where " dual AGN" were decided a more likely scenari o 
than outflows or rotating disks (IKomossa et al.ll2008l ; IComerford et al.ll2009l ; iLiu et al.ll2010l ). 



Crenshaw et al.l ( I2010af ) showed that velocity offsets between the narrow emission lines 



and host galaxy lines can be explained by the combination of mass outflows and dust ex- 
tinction in AGN, or even an asymmetric distribution of outflowing gas. They also sug- 
gested that these explanations could possibly work for double-lined AGN, such as Mrk 78, 
as well. Previous studies had already suggested that radial motion via inflo ws or outflows 



could be responsible for the narrow line asymmetries and velocity shifts ( jWhittld 11985 
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Dahari fc De Robertia Il988l ; iDe Robertis Shawl Il990l ). including long-slit ground-based 



observations that gave evidence for the production of asymme t ries via outflows and dust 



obscuration (IStorchi-Bergmann et al 



Ramos Almeida et al 



2006 



1992 



Arribas et a 



19961 : IChristopoulou et all Il997l : 



Rodrfguez-Ardila et al.ll2006l ). In this paper, we use high reso- 



lution spectroscopy and imaging to show that the velocity offset seen in Mrk 78 can indeed 
be explained by mass outflows in the NLR. 



2. Observations 

We obtained archival Space Telescope Imaging Spectrograph (STIS) long-slit spectra 
and Faint Object Camera (FOC) images of Mrk 78 from the Multimission Archive at the 
Space Telescope Science Institute (MAST). Four sets of spectra were taken between 1998 
February 28 and 1998 March 01 under Hubble Program ID 7404 (M. Whittle, PI) with a 
52" x 0'/2 slit using a medium-dispersion G430M grating (4950 - 5240A), which includes 
[O III] A5007. The spectral resolution was 0.56 A(~ 30 km s" 1 FWHM), with an angular 
resolution of 0'.'051 per pixel in the cross dispersion direction. Of the four observations, three 
were taken at a position angle of 88° and the fourth was taken at a position angle of 61°. 
Additional details of each observation are given in Tables 1 and 2. 

Figure 1 shows an extended continuum STIS acquisition image of Mrk 78 with addi- 
tional, contoured [O III] and UV FOC images and positions for the four STIS slits illustrated 
on the [O III] image. The slits are offset from the optical continuum peak, signified as a 
"+", which was located by alig ning the [O III] and continuum FOC images taken on the 



same date (Schmitt et al. 



Ramos Almeida et al.l ( 120061 ). 



20031) and is also coincident with the infrared peak as observed by 



The STIS spectra were processed using IDL software developed at NASA's Goddard 
Space Flight Center for the STIS Instrument Definition Team. Cosmic-ray hits were removed 
before further processing. The zero point of the wavelength scales were corrected using 
wavelength-calibration exposures taken after each observation. Finally, the spectra were 
combined, geometrically rectified, and flux calibrated to provide spectral images that have a 
constant wavelength along each column and display fluxes in units of erg s _1 cm~ 2 A -1 per 
cross- dispersion pixe l. The calibrated spectral images for each slit position can be found in 



Whittle et al.l (120021 ). 



We used the STIS spectra to compare our results with those taken from previous ground- 
based spectra, obtained through apertures that included most, i f not all, of the NLR in 



Mrk 78 at moderate spectral resolution (80-230 km s" 1 FWHM) (IVrtilek fc Carletonlll985 
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Nelson fc Whittle! Il995h . The STIS datasets consist of hundreds of spectra along each slit, 
allowing us to resolve individual emission- line knots that are responsible for a majority 
of the emission seen in NLRs. To simulate the ground-based spectra, we combined all 
of the STIS spectra along each slit position and then combined the spectr a from each slit 



excludi ng slit D as it overlaps the other three slits, similar to the procedure in lCrenshaw et al. 



( 12010bl ). Figure 2 shows the final co-added spectrum in terms of flux as a function of the 
radial velocity of the [O III] line, with respect to the systemic velocity calculated from our 
redshift. The significantly separate d [O HI] A5007 pea k s agre e with those from the moderate- 
resolution ground-based spectra of iNelson &: Whittle! (119951 ). with the higher peak being at 
near systemic velocity (v r = —50 km s _1 ) and the secondary red peak at v r = 800 km 
s _1 , meaning any [O III] emission not covered by the STIS long-slit observations does not 
significantly alter the profile shape or structure. 

To determine the kinematics of the NLR, we incorporated the processed spectral images 
into a program that fits [O III] A5007 emission-line components with G aussians oyer a n 
average continuum, taken from available line- free regions in the spectra (IDas et al.l 120051 ). 
In order to get the maximum number of measured lines per data set, spectra were smoothed 
with a boxcar average using a width of 7 pixels. This smoothing allowed for the fitting of 
extended, low flux [O III] emission lines, without inhibiting us from deblending multiple 
[O III] lines dependent on different kinematic components along the slit. Calculating the 
central peak location for each Gaussian allowed us to determine the central wavelength 
for the corresponding [O III] line, which in turn was used to calculate its Doppler shifted 
vel ocity. There wer e three sources of uncertainty in our velocity measurements, as detailed 
by lDas et al.l (120051 ). The first is that the measured emission lines are not perfect Gaussians, 
the second comes from emission cloud displacements from the center of the 0'.'2 slit in the 
dispersion direction, and the third comes from noisy spectra. Errors were converted to 
velocities and added in quadrature, which produced a total maximum error of ± 60 km s _1 . 
Extremely noisy spectra (spectra without detectable emission > 3a) at locations > 3" — 3'.'5 
to the west and > 3'/75 (4'.'25 for slit D) to the east of the nucleus, were not fitted. Subtracting 
the systemic velocity of the galaxy (cz = 11037 km s _1 ) from velocities calculated via [O III] 
line centroid shifts, we were able to determine the radial velocities along each slit in the rest 
frame of Mrk 78. 



3. Results 



Figures 3 and 4 show the rest frame radial velocities, FWHM, and normalized fluxes 
for the measured [O III] A5007 emission lines from the four slit positions. Kinematics from 
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the lower spectral resolu ti on IR spectrum of iRamos Almeida et al.l (120061 ) and the spectral 
images of IWhittle et al.l ( 120021 ) agree with our measurements. In general, the structure 
of the emitting region is clumpy in nature, as seen in the large number of flux peaks in 
each slit that can be matched to individual bright knots in the [O III] image of Figure 
2. While each knot appears to have an individual peculiar velocity, the velocity field as a 
whole is relatively well organized. High radial velocities on both sides of the nucleus show 
the distinct asymmetrical red/blue shifts characteristic of biconical outflows seen in several 



other l ow-redshift Seyfert galaxies (IRuiz et al.l l2001t iDas et al.ll2005l 120061 ; ICrenshaw et al. 



2010al ; iFischer et al.ll2010l ). Emission eastward of the optical peak is dominantly redshifted 



from systemic velocity (henceforth 'red lobe') while emission westward of the optical peak 
is dominantly blueshifted (henceforth 'blue lobe'). Breaking down our simulated ground 
based spectrum of Figure 2 into the individual slits A, B, and C, we find that the redshifted 
velocity peak is primarily due to a specific high-flux knot, best sampled in slit B. Figure 1 
contains a circle around the location of the redshifted knot and Figure 3 has a bracket in slit 
B, located over the high-flux region which depicts the knot. 



Models 



To generate kinematic models, we assume the ionizing radiation responsible for the NLR 
is biconical in nature, as it is the simplest geo metric shape produced by a central obscuring 



torus. As with our previous kinematic models (IDas et al. 



2005, 2006; Crenshaw et al. 2010a 



Fischer et al.l 120101 ) . we also initially assume that the bicone is hollow, where the separate 
NLR emission co mponents represe nt opposite sides of the bicone. We employed our kinemat- 
ics program from ( IDas et al.ll2005l ). which allows us to recreate the observed radial velocities 
along a fixed slit position by altering various parameters of our model bicone outflow^ Our 
models include seven alterable parameters with initial values taken from the [O III] image in 
Figure 1 (deprojected height of bicone [2^], outer half-opening angle [S m(B ], and position 
angle) and the kinematics from Figures 3 and 4 (maximum velocity [iWr] and turnover ra- 
dius [r t ]), leaving the inner half-opening angle (8 min ) and inclination of the bicone axis out of 
the plane of the sky to be altera ble parameters free of restriction. We adopted our previous 
velocity law (IFischer et al.ll2010l ). a linear increase starting with zero km s _1 at the nucleus 



1 Previous studies of Mrk 78 have suggested th at the NLR outflow is being radially accelerated via radio 
jets |Whittle fc Wilson! booi IWhittle et al.ll2005h . While our model does not include any specific source of 
acceleration, instead simply simulating general outward acceleration, we do not detect considerable differ- 
ences in the velocities o r dispe rsions of emission-line knots near radio knots from others in the overall outflow 



noted in IWhittle et al] (|2002f ). 
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and subsequent linear decrease at r t ending at zero km s 1 , because it is the simplest law 
that matches the observations. 

The FOC [O III] image of Figure 1 shows that the NLR is not symmetric on either 
side of the nucleus. The redshifted, NE cone is less dispersed, more conical in nature, and 
features cleaner kinematics from which we can discern separate kinematic components due to 
outflow emission in the fore and back side of the cone. Conversely, the blueshifted, SW cone 
appears to have a larger opening angle, with floculent, knotty emission and resulting intricate 
kinematics. Though using initial parameters solely from the NE cone was enticing, we found 
that the bicone did not envelop the entirety of the west lobe. Thus, we expanded the model 
bicone outer opening angle to envelop all NLR emission in the [O III] image in Figure 2. We 
chose Omax — 30° as our initial opening angle and used the kinematic components of the NE 
cone as the basis for the majority of our fitting. The initial inclination parameter relied on the 
kinematics seen in Figures 3 and 4, as either side was either mostly redshifted or blueshifted, 
and the imaging in Figure 1 as the extended biconical structure on either side of the nucleu s 
would be less pronounced if we were to be looking into the bicone (ISchmitt fc Kinneylll996l ). 
With these guidelines, we determined an intermediate angle was most appropriate for a 
starting position. To find the most successful model, we altered parameters across all data 
sets while applying specific centroid offsets and position angles for each slit until a specific set 
of free parameters created an acceptable fit to all four data sets. Final values for all model 
parameters are provided in Table 3. Interestingly, the final opening angle and position angle 
are nearly identical to the apparent, initial parameters taken fro m Figure 1 (9 mnT = 30° and 
P. A. = 70° respectively). In co mparison to our previous studies ( IRuiz et al.ll2001t iDas et al. 



20051 . 120061 ; iFischer et al.ll2010[ ). Mrk 78 has an extremely large NLR as well as a narrower, 
thicker (greater difference between 6 max and m in) bicone than all other targets thus far. 

Figure 5 displays a good qualitative fit between the model created by our final set 
of parameters and the data. The model matches fairly well to the east of the nucleus as 
kinematics for both sides of the eastern lobe are discernible in multiple slits. One discrepancy 
worth discussing is the lack of model data (shaded regions) mainly to the west of the nucleus 
in Slits B and C between — 0'/25 to 0'.'5 and — 0'.'5 to 1'.'5 respectively. As our kinematics 
model creates a well defined geometrical model bicone with a sharp apex, incorporating 
slit positions that are transversely offset from the nucleus will generate models that lack 
data near the nucleus as the slit does not pass over the model bicone. This is due to the 
simplicity of our model as we assume that the bicone has a sharp apex, which is obviously 



not t he case in this and other Seyferts (lAntonuccilll993l ; ICrenshaw et al.ll2010bl ; IFischer et al. 



20101 ) . Thus, slits B and C lack modeled kinematics in these regions as the actual emission 
lies outside our model bicone geometry. 
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5. Discussion 

The available evidence indicates that Mrk 78 is not a double AGN. Observing the 
variance in flux across the slit, we see many bumps and peaks that are due to large, bright 
knots of emission-line gas, each with their own peculiar velocity inside the general outflow 
pattern. It is one of these knots, a bright red-shifted cloud traveling at ~ 800 km s _1 inside of 
the ionizing bicone, that is the source of the offset radial velocity pe ak. The cloud appears to 



fit well in our modeled outflow pattern, supporting the suggestion by lCrenshaw et al.l (J2010bJ) 



that uneven outflow can account for such dual-peaked profiles, and thus there is no need to 



the NLR in Mrk 78 actuallv look rather similar to those in NGC 4151 (Hutchines et al. 19981; 


Crenshaw et al. 2000: Das et al. 20 


0, 


5), NGC 1068 


Crenshaw & Kraemer 


20001: feas et al. 


2006. 


2007). Mrk 3 


Ruiz et al. 2001 


: ( 


Drenshaw et al. 


2010a), and Mrk 573 ( 


Schlesineer et al. 


2009: 


Fischer et al. 


2010). For each of these tareets. including Mrk 78. we have successfully 



employed our biconical outflow models to find a good match to the overall flow patterns as a 
function of distance from the central AGN. Though there are some inconsistencies between 
the model and the data that are likely due to the simplified geometry of our bicone model 
and pecular velocities of individual emission-line knots, radial outflow of [O III] clouds can 
explain the general, large-scale flow pattern of the NLR in Mrk 78, as well as the double 
peaked emission. 

Does obscuration fro m dust contribute to the creation of the dual-peaked radial veloc- 
ity profile as mentioned in ICrenshaw et al.l ( )2010bl )? Ground based images of Mrk 78 in B 
and I give a position angle and inclination o f the host galaxy as 84° and 64° respectively 
( jSchmitt fc Kinneyl I200CH : Kinney et al.l 120001 ) . Adding this geometry to a geometric bicone 
model, we can test to see if a smooth distribution of dust contributed by the disk could be 
responsible for the radial velocity profile. Figures 6 and 7 show simple, three-dimensional 
geometric models of the ionizing bicone's outer surface using 9 max , position angle, and in- 
clination values from our final kinematics model and an additional two-dimensional galactic 
disk using the above position angle and inclination values. Though it is unknown which side 
of the disk is closer, either geometry in combination with our final bicone model results in a 
radial velocity profile that is either symmetric or asymmetric towards the blue. The latter 
case can be seen in Figure 6 as the redshifted cone should be more extincted by dust in the 
galactic plane. This scenario does not fit our data as the blue-shifted, likely less obscured 
lobe would then contribute a majority of the total emission-line flux. Thus, it is unlikely 
that a dusty disk would be responsible for the given profile. 

There is a prominent dust lane that sweeps over the center of the AGN. Each image 
in Figure 1 shows a distinct band from southeast to northwest across the nucleus with a 
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position angle of ~123° southeast to northwest that has flux much lower than the emission 
surrounding it. More quantitatively, we obtained an additional STIS G750M (5450 - 10140A) 
spectrum from MAST taken along the same position as the G430M slit D in order to measure 
the flux ratio of Ha/Hf3 over the nucleus, and calculated a reddening of E(B — V) — 
0.78 (assuming a galactic reddening curve and intrinsic HajB.fi = 2.9). The reddening is 
significant enough to extinguish portions of the inner NLR, possibly suppressing high- velocity 
dispersion emission from the bicone at ~ km s _1 and making the double line profile more 
distinct. However, this reddening is not enough to extinguish the AGN broad-line region 
(BLR). As such, it is likely that this dust lane must be additional obscuration external to 
the "torus". 

Further evidence for a dust lane can be seen in the blue-shifted southwest lobe from 
the imaging shown in Figure 1. From the end of the absorbing dust lane, an arc of emission 
is visible along the north rim of the lobe ~ 2" from the "+" marking the continuum peak 
location. It remains unclear whether the dust lane crossing across the nucleus arcs over the 
southwest lobe and continues to absorb emission or ins tead arcs into the ionizing bicone and 



becomes illuminated, similar to Mrk 3 and Mrk 573 ( Crenshaw et al.ll2010at iFischer et al. 



20101 ). However, both situations would help explain the asymmetrical, patchy emission seen 
in the blue-shifted emission lobe. Figure 6 depicts a galactic disk with the north side of 
the disk coming out of the page, which allows for a potential interaction between the disk 
and ionizing bicone. A full understanding of this extended morphology will require more 
thorough imaging before a more accurate claim can be made. 



6. Conclusions 

We have analyzed STIS long-slit G430M spectra of the inner emission knots and nucleus 
of the Seyfert 2 galaxy Mrk 78. We generated kinematic and geometrical models of the NLR 
and ENLR that successfully fit the data, with the final parameters of Mrk 78's outflow bicone 
given in Table 3. While it is not clear at this time how the blue-shifted lobe is interacting 
with the dust lane lying over the nucleus, it is likely that that interaction can explain the 
dispersed morphology present in the given FOC and STIS imaging. 

We find that the double peaked emission lines seen in ground-based spectra of Mrk 
78 are due to an asymmetric distribution of outflowing gas in the NLR. Specifically, the 
red peak at +800 km s _1 is due to a bright emission-line knot that is redshifted with 
respect to the systemic velocity of the host galaxy. The bright, redshifted knot follows the 
same flow pattern as the other knots, and thus, there is no need to assume that a second 
AGN resides within Mrk 78. Furthermore, all Seyferts that we have modeled to date show 
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outflows (IRuiz et al.l 1200 ll ; iDas et al.l 12005k l2006t iFischer et al.l 120101 ) , and thus asymmetric 
distributions of outflowing gas and/or patchy dust extinction maybe a general explanation 
for double-peaked emission-line profiles. With nothing unusual separating Mrk 78 from other 
AGN, it is possible that many double-peaked profiles may be due to outflows in a single AGN 
as shown by our models. 
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Table 1. HST/ST1S G430M Observations of Mrk 78 



Slit 


Exposure 
(s) 


P.A. 

(degrees) 


Offset" 
(arcsec) 


A 


1727 


88 


0.125 


B 


1938 


88 


-0.27 


C 


2052 


88 


-0.55 


D 


1800 


61 


-0.05 



a Cross-dispersion offset from the op- 
tical continuum peak. A positive value 
corresponds to a displacement south of 
the slit P.A. axis in Figure 1. 



Table 2. i/ST/FOC/STIS Imaging Observations of Mrk 78 



Instrument 


Filter 


Exposure 
(s) 


STIS 


MIRVIS 


120. 


FOC 


372M 


896. 


FOC 


550M 


1196. 


FOC 


502M 


801. 
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Table 3. Best Fit Model Parameters for Mrk 78 a 



Parameter Values 



P.A. 65° 

i 30°(SW) 

max uu 

drain 10 

v max 1200 km s" 

^max 

3200 pc 

r t 700 pc 



a The letters in parenthe- 
ses indicates the side closest 
to us. 
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Fig. 1. — Composite of HST STIS acq uisition continuum i mage (left) with additional en- 



larged [O III] (F502M filter, as seen in ISchmitt fc Kinnevl (119961 )) and UV (F372M filter, 



mainly [Oil] emission) FOC images (middle and right respectively). The dust lane sweep- 
ing over the AGN is located at a position angle of ~123° southeast to northwest. STIS 
slit positions have been superimposed and the red-shifted emission knot responsible for the 
secondary radial velocity peak is circled on the [O III] image. 

Fig. 2. — Summation of spectral flux from slits A, B, and C as a function of radial velocity 
for the [O III] A5007 emission line with respect to the systemic velocity. Double peaked 
[O III] A4959 emission-line profile is seen to the left. 

Fig. 3. — Radial velocities (top), FWHM (middle), and normalized total flux (bottom) of 
[O III] measurements for slits A and B. The bracket in slit B shows the location of the high 
flux cloud responsible for the secondary radial velocity peak. 

Fig. 4. — Radial velocities (top), FWHM (middle), and normalized total flux (bottom) of 
[O III] measurements for slits C and D. 

Fig. 5. — Shaded regions represent the kinematic model chosen as the best fit for our radial 
velocity data set for all four slits. Parameters used to create this model are given in Table 
1. Lighter regions represent kinematics for the side of each cone near parallel to our line of 
sight, darker regions represent kinematics for the side near perpendicular to our line of sight. 

Fig. 6. — Geometric model of the NLR based on parameters from Tab le 1, shown as viewed 
from Earth. Inner disk geometry is taken from iKinney et al.l (120001 ) depicting the top of 
the disk coming out of the sky. The yellow bar represents the position angle of the dust 
lane traveling over the nucleus. The bright redshifted knot is in the northeast portion of the 
bicone. 



Fig. 7. — Same as Figure 6, except the inner disk geometry depicts the bottom of the disk 
coming out of the sky. The yellow bar represents the position angle of the dust lane traveling 
over the nucleus. 
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Fig. 5. 
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Fig. 7. 



